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Abstract 

i-C ■ If the positron fraction and combined electron-positron flux excesses recently ob- 

served by PAMELA, Fermi and HESS are due to dark matter annihilation into lepton- 
rich final states, the accompanying final state radiation (FSR) photons may be detected 
by ground-based atmospheric Cherenkov telescopes (ACTs). Satellite dwarf galaxies 
in the vicinity of the Milky Way are particularly promising targets for this search. 
We find that current and near-future ACTs have an excellent potential for discover- 
ing the FSR photons from dwarfs, although a discovery cannot be guaranteed due to 
large uncertainties in the fluxes resulting from lack of precise knowledge of dark matter 
distribution within the dwarfs. We also investigate the possibility of discriminating be- 
tween different dark matter models based on the measured FSR photon spectrum. For 
typical parameters, we find that the ACTs can reliably distinguish models predicting 
dark matter annihilation into two-lepton final states from those favoring four-lepton 
final states (as in, for example, "axion portal" models). In addition, we find that the 
dark matter particle mass can also be determined from the FSR spectrum. 



1 Introduction 

Recent measurements of the positron fraction in cosmic rays in the 10-80 GeV range by 
PAMELA [1] and the combined electron-positron flux up to the TeV scale by Fermi [2] and 
HESS [3] indicate excesses inconsistent with conventional astrophysical background from 
cosmic ray nuclei interactions with the interstellar medium. Although the possibility that the 
origin of these excesses could lie within conventional astrophysics, such as nearby pulsars [4] , 
remains, there is widespread optimism that these excesses could be the first indirect evidence 
of dark matter annihilation in the galaxy. 

The simple hypothesis that dark matter consists of weakly interacting massive particles 
(WIMPs) yields the correct current dark matter density, and is supported by many popular 
extensions of the standard model of particle physics. It can also fairly naturally explain these 
excesses. The WIMP framework requires the dark matter particles to have a (thermally 
averaged) annihilation cross-section of (av) ~ 3 x 10 -26 cm 3 s _1 at freeze-out, and a mass in 
the range from 10 GeV to 10 TeV. PAMELA, Fermi, and HESS excesses can be explained, in a 
way consistent with all other astrophysical observations, by dark matter with an annihilation 
cross section of (av) ~ 3 x 10~ 23 cm 3 s _1 in the Milky Way, and a mass of 1-3 TeV, provided 
annihilation is predominantly into electrons or muons [5, 6]. 1 The (9(100-1000) discrepancy 
in annihilation cross section can be bridged by incorporating Sommerfeld enhancement [9, 10] 
and/or a "boost factor" associated with small-scale dark matter clumps in the Milky Way 
halo 2 [11]. 

If the observed excesses indeed have a dark matter origin, accompanying signals are ex- 
pected in the form of energetic gamma rays. For leptophilic dark matter, these are produced 
in two ways: inverse Compton scattering (ICS) of starlight and CMB photons off energetic 
leptons from dark matter annihilation, and final state radiation (FSR) from the annihilation 
process [12, 13, 14, 15, 5, 16]. We will focus on the FSR photons in this paper. The FSR 
component is dominant at energies close to the dark matter mass. It is also independent 
of most of the astrophysical uncertainties that plague the calculation of positron and ICS 
photon fluxes. Finally, if WIMP annihilation primarily results in two-body leptonic final 
states, the associated FSR spectrum has a sharp "edge" feature at the energy equal to the 
WIMP mass, enhancing its visibility over astrophysical backgrounds and allowing for precise 
mass measurement if an edge is observed [17]. 

Dark matter annihilation rate, and with it the FSR photon flux, are proportional to 
the square of the dark matter density, and regions with enhanced density should provide 
the strongest flux. The largest concentration of dark matter in our cosmic neighborhood 
is of course in the central region of the Milky Way. However, fluxes of hard gamma rays 
from astrophysical sources located near the galactic center are large and poorly understood, 



lr The most significant constraint on these models currently arises from the first-year Fermi-LAT diffuse 
gamma ray map [7]. An isothermal-like dark matter density profile in the Milky Way is required to satisfy 
this constraint. The constraint is especially important for the muon final states [8]. 

2 The degree of clumps that would provide all of the required enhancement, however, is disfavored by 
numerical simulations [111. 



providing a serious background to dark matter searches in that region. In contrast, dwarf 
galaxy satellites of the Milky Way, which are highly dark matter dominated, are largely free 
of astrophysical backgrounds, and many of them lie away from the galactic center [18]. This 
led many authors to advocate dwarfs as promising sources for gamma ray signals of dark 
matter [19, 20, 16]. 

The Fermi Large Area Telescope (LAT) can cover the whole sky continuously, detect 
photons down to 100 MeV energy, has superior energy and angular resolution, and is free of 
atmospheric background. However, as dark matter gamma ray signals from dwarf galaxies 
are mainly constrained by low statistics, atmospheric Cherenkov telescopes (ACTs) offer a 
distinct advantage: the typical ACT effective area is O(10 9 ) cm 2 , versus the ~ 10 4 cm 2 
effective area of Fermi, offsetting the disadvantages of limited observation time, poorer res- 
olution, and atmospheric background. 3 An additional shortcoming of the ACTs in many 
searches is their high energy threshold, typically around 200 GeV. However, as the FSR pho- 
ton spectrum extends all the way to the energy corresponding to the WIMP mass, which is 
required to be in the 1-3 TeV range to explain PAMELA, Fermi and HESS electron/positron 
anomalies, this limitation does not play a crucial role in our case. 

Observation of dwarf galaxies with ACTs for possible signals of dark matter annihilation 
has been carried out on numerous occasions, without any positive results so far. These 
searches already put some bounds on leptophilic models of dark matter, as we discuss in 
Section 3. (Some of these bounds have also been considered in Refs. [14, 16, 22].) However, 
the bounds depend sensitively on the distribution of dark matter inside the dwarfs, which is 
at present poorly known. As a result, none of the models we examine is currently ruled out. 
Future searches can benefit from recent discoveries of new dwarf galaxies by the Sloan Digital 
Sky Survey [23] - in particular, the dwarf galaxy Segue 1 [24] - as well as rapid advances 
in the field of ACTs - MAGIC II recently became operational [25], while the Cherenkov 
Telsecope Array (CTA) [26] and the Advanced Gamma-ray Imaging System (AGIS) [27, 28], 
which will provide an order of magnitude improvement over current ACTs, are due to be 
completed in the next few years. We discuss the prospects of detection in Section 4. We 
find that the ongoing and near-future ACTs will cover most of the flux range predicted 
by the leptophilic models, although unfortunately the large uncertainties on dark matter 
distribution in the dwarfs do not allow us to conclude that the discovery is guaranteed. 

Inspired by the discovery prospects, in Section 5 we investigate the important question 
of whether the ACTs can measure the FSR photon spectrum precisely enough to distin- 
guish between different dark matter annihilation channels predicted by different models of 
dark matter. We find that, for typical dark matter and astrophysical parameters, model 
discrimination prospects are good: in particular, the ACTs can reliably distinguish models 
which predict two-body final states from those favoring four-body final states. The latter 
possibility is common in models where annihilation proceeds via a scalar "portal" particle; 
the same particle may give rise to the Sommerfeld enhancement of the cross section [10, 29]. 
Discrimination between the different four-body final states (4e vs. 4/j,) is more challenging, 



3 Fermi LAT may be used to constrain properties of dark matter with gamma rays originating from the 
galactic center, see Ref. [21]. 



since the FSR spectra in these two cases are quite similar. In addition, we find that spectral 
fits provide accurate determination of the dark matter mass and annihilation cross section. 

2 The Final State Radiation Spectrum 

Direct annihilation of WIMPs into photon pairs cannot occur at tree level, and thus has 
strongly suppressed cross section. In contrast, final state radiation (FSR) of photons can 
occur whenever WIMPs annihilate into electrically charged particles, i.e. in the vast majority 
of possible SM final states (the only exception being all-neutrino channels). Most of the 
energetic FSR photons are approximately collinear with the emitting charged particle. In 
this kinematic regime, factorization theorems ensure that the energy spectrum of the FSR 
photons is to leading order independent of the details of the WIMP annihilation process, so 
that quasi-model-independent predictions can be made. 4 Moreover, for WIMP annihilations 
into two-body fermionic final states, e.g. XX ~~ *" e + e~, the FSR spectrum has a characteristic 
"edge" feature at the energy equal to the WIMP mass, providing a powerful signature to 
discriminate it against the poorly known astrophysical background. All these features make 
FSR photons a promising candidate for dark matter searches, as emphasized in Ref. [17]. 

Specifying the final state of the WIMP annihilation process completely fixes the FSR 
spectrum. In this paper, we will focus on the final states containing charged leptons. We 
will assume that there is a single WIMP particle, x-, with the mass in the 1-3 TeV range, 
that can annihilate into leptons either directly, or via an intermediate "portal" particle <p, 
with mass of order 1 GeV. Specifically, we consider the following three annihilation channels: 

1. Model A: XX ~> ^ + ^ ■ 

2. Model B: XX ->• 00 ->■ 4e. 

3. Model C: XX ->■ 00 ->■ V 

These parameters and final states are favored by the combination of PAMELA, Fermi and 
HESS data with bounds on gamma ray flux from the galactic center and other astrophysical 
constraints [5, 6]. 

For dark matter annihilating into a lepton-antilepton pair 11 as in model A, the FSR flux 
(within the leading- log approximation) is [17] 



d ®FSR ^ /(<rt/)\ /lOOGeV\^ . , (Am 2 Jl-x 



where x = 2E 1 / \fs = E 7 /m x , $o = 1-4 x 10 14 cm 2 s 1 GeV 1 , F(x) is the splitting function 
for fermions 

F(x) = i±(iz£>!, (2) 



X 



4 This is analogous to the quasi-model-indcpendcnt predictions for WIMP production at colliders in as- 
sociation with a photon or a gluon [30]. 



and J is the dimensionless astrophysical factor 5 

J= — \ — [ *] L, L=[dn[ p 2 dl. (3) 

8.5 kpc V0.3 GeV/crn 3 / J Ji.o.s. w 

For XX — >• 00 — >■ 4/, as in models B and C, the FSR spectrum is [17] 

d®FSR ^ f(<rv)\ /l00GeV\ 3 2 2-:r + 2:rlog:r-:r 2 ^ /ro|\ J 



<ix y lpb J \ m x J x \ raj j 

For the remainder of this paper, we use (av) = 3 x 10~ 23 cm 3 s _1 (corresponding to, for 
example, an unboosted do = 1 pb and a factor of 1000 from Sommerfeld enhancement), 
m x = 3 TeV, and m^ = 1 GeV 6 unless stated otherwise. These parameters provide consistent 
fits to PAMELA, Fermi and HESS data for all three models at hand [5, 6]. While lower 
velocities in the dwarf galaxies can lead to greater Sommerfeld enhancement, additional 
constraints from CMB observations [32] suggest that they cannot be much larger than the 
value we have chosen at this mass scale. Note that the use of leading-log approximation 
assumed in Eqs. (??) is well justified for these parameters: the largest correction beyond the 
leading-log occurs in model C, and is of order O^-^-) 2 ) ~ 0.04. 

For models A and C, annihilation is to muons, and final state radiation form the subse- 
quent decay of the muon also needs to be taken into account. The FSR photon spectrum 
from radiative muon decay in the muon rest frame is [33, 15]: 

d% al-w{ /n n , o n , N1 1-w 17 23w 101w 2 55w 3 \ 

— 1 = 3 - 2w + Aw 2 - 2w 3 ) In — + — — - + — - , 5 

aw Sit w \ r 2 6 12 12 J 

where r = m 2 e /m 2 , w = 2E 1 /m^. (An approximation r C 1 has been made in this formula.) 
If the muons come directly from dark matter annihilation, this spectrum needs to be boosted 
to the dark matter rest frame: 

^I = 2 /'feI^, (6) 

ay J y w aw 



5 Eq. (1) and similar formulas below are equally applicable for decaying dark matter, with an appropriate 
redefinition of the J factor. However, we do not consider decaying dark matter in this paper since the 
resulting FSR signals from dwarf galaxies are generally too weak to be detected. See Ref. [31] for a discussion 
of discriminating between decaying and annihilating dark matter based on gamma ray observations. 

6 A 1 GeV portal particle is only consistent with observations if it does not have significant couplings 
to quarks, since otherwise <p decaying to neutral pions produce excessive gamma-ray flux from the galactic 
center. We assume such a scenario in this paper. An alternative is to consider a portal particle that does 
couple to quarks (such as, for example, an extra U(l) gauge boson mixed with the photon), but is too light 
to decay to pions. In this case, decays to muons are only allowed if m^ is in a narrow range about 220-250 
GeV, which seems unnatural. Let us note, however, that the formalism of this section would apply to this 
case as well, with a possible exception of model C where (m fM /m ( p)' 2 suppressed corrections would need to be 
included in the FSR photon spectrum. 



where y = E 1 /m x . When annihilation is through an intermediate particle, Eq. (5) needs 
to be boosted twice to get to the dark matter rest frame [16]: first to rest frame, 

- 2 / dw — —*, (7) 



dx J2x/(i+(3) w dw 

where x = E^/m^, followed by a boost to the x res t frame, 

d$ r 1 1 dN 

^2 = 2 / dx-^p, (8) 

dy J y x dx 

where y = E 1 /m x as before. The overall FSR signal is obtained by adding this contribution 
to FSR off the muons. Typically, FSR off muons from the annihilation process remains 
dominant unless m^ ~ m M [16]. 

3 Dwarf Galaxies and ACTs 

In this section, we will summarize the relevant properties of the known dwarf galaxies that 
are potential sources of observable FSR photons, as well as the basic parameters of existing 
and near-future atmospheric Cherenkov telescopes (ACTs). We will also discuss the bounds 
placed on the models at hand by the previous ACT searches for anomalous gamma rays from 
dwarfs. 

3.1 Dwarf Galaxies 

Dwarf satellite galaxies, made up almost entirely of dark matter, provide excellent prospects 
for detecting dark matter signals. After the galactic center - where astrophysical backgrounds 
are large and poorly understood - these dwarfs have the largest dark matter density in our 
galactic neighborhood. In addition, they have no detected neutral or ionized gas, minimal 
dust, no magnetic fields, and little or no recent star formation activity [18]: in other words, 
one should expect essentially no astrophysical background associated with the dwarfs. 

Table 1 lists the major known dwarf galaxies that are promising candidates for dark 
matter detection, with their astrophysical factors as defined in Eq. (3) and their galactic 
coordinates. The astrophysical factors are taken from Ref. [16], where they are calculated 
for AQ ?» 2.2 x 10~ 5 sr, which corresponds to the angular resolution of the HESS telescope, 
A6 = 0.15°. These dwarfs have already been observed by one or more ACTs, except Segue 1, 
which is a recently discovered candidate [24] that seems to be particularly promising because 
of its large dark to luminous matter ratio, relative proximity, and high altitude. Segue 1 has 
been observed by Fermi [34] with null results. We will not use Sagittarius in our analysis as it 
is being tidally disrupted by the Milky Way, making the estimates of the L factor unreliable. 
In addition, it also suffers from a larger background since its direction is close to that of the 
galactic center. 

The Sloan Digital Sky Survey [23] has recently discovered many new dwarf galaxies. Since 
only a small region of the galactic neighborhood has been completely surveyed, it seems likely 



Dwarf 


log 10 (Lx GeV- 2 cm 5 ) [16] 


galactic coordinates (/, b) [35] 


Sagittarius 


19.35 ±1.66 


5.6°, - 14.2° 


Draco 


18.63 ±0.60 


86.4°, 34.7° 


Ursa Minor 


18.79 ±1.26 


105.0°, 44.8° 


Willman 1 


19.55 ±0.98 


158.6°, 56.8° 


Segue 1 


20.17 ±1.44 


220.5°, 50.4° 



Table 1: Dwarf galaxy candidates for dark matter indirect detection. 



that several hundred more low-luminosity, dark matter dominated dwarf galaxies might still 
be discovered in the future [36]. It would be straightforward to apply the analysis of this 
paper to any promising new dwarf that may be discovered, once the distribution of dark 
matter is mapped out to allow for at least an approximate determination of its L factor. 

3.2 Atmospheric Cherenkov Telescopes 



ACT 


Effective area (cm 2 ) 


energy resolution e 


MAGIC [37] 


10 9 


0.20 


HESS [38] 


7x 10 8 


0.15 


VERITAS [39] [40] 


7x 10 8 


0.15 


CANGAROO III [41] 


10 9 


0.15-0.30 


STACEE [42] 


10 8 


0.25 


MAGIC-II [25] [43] 


10 9 


0.10-0.15 


AGIS [27] [28] 


10 10 


? 


CTA [26] 


10 10 


0.10 



Table 2: Some current and near-future ACT's relevant for dwarf galaxy observations. 
MAGIC II started taking data in 2009. CTA is expected to be operational by 2015. Con- 
struction of AGIS could start between 2011-2013. 



Since gamma ray fluxes from FSR in dark matter annihilation in dwarf galaxies are rather 
weak, the key parameter governing the sensitivity of a given telescope in this search is its 
effective area A e g. In addition, an important ACT parameter is the energy resolution of the 
instrument e. The probability to assign an energy E to a primary gamma ray of true energy 
E' can be approximated with a Gaussian as [20] 



R e (E - E') 



1 



2nE't 



cxp 



(E-E 
2e 2 E' 2 



A2 N 



(9) 



This parameter is particularly important for measuring the gamma ray spectrum, which 
can be used to identify the dark matter model if an FSR signal is observed (see Section 5). 
Table 2 lists current and future ACT's relevant for indirect detection of dark matter. The 
listed effective area values are averages in the several hundred GeV to 1 TeV range; more 
detailed plots of effective area as a function of energy for individual instruments can be 
obtained from the respective references. Instead of using telescope-specific, energy dependent 
effective areas, we ignore energy dependence and use typical parameters in this paper so 
that our estimates remain more generally applicable. Unless stated otherwise, we assume an 
angular resolution of A6 = 0.15°, corresponding to a solid angle coverage AQ « 2.2 x 10~ 5 
sr. The typical instrumental energy threshold for ACT's is about 200 GeV; however, in 
some cases imposing a higher threshold in the analysis can improve the sensitivity of FSR 
searches, as we discuss below. 

3.3 Previous Observations and Upper Bounds 

Previous ACT observations of dwarf galaxies are listed in Table 3. No significant signals 
were observed, resulting in the listed upper bounds on high energy gamma ray flux from 
dark matter annihilation in these galaxies. The VERITAS upper limits for the three dwarfs 
were at the level of 1% of the Crab nebula [44], which is the value listed in the table. 
As mentioned earlier, we will not use the Sagittarius bound in our analysis; however, it is 
worth mentioning that earlier studies have found dark matter interpretations of PAMELA 
and Fermi anomalies to be incompatible with this bound for optimistic halo profiles of 
Sagittarius [14]. Fermi LAT telescope has also observed several dwarf galaxies, placing 
flux upper limits of 0(1O~ 10 ) cm _2 s _1 above 1 GeV and 0(1O~ 9 ) cm _2 s _1 above 100 MeV; 
individual bounds are listed in [45]. 



Dwarf 


observed by 


total exposure 
time (hrs) 


upper bound on photon flux (cm 2 s r ) 


Sagittarius 


HESS 


11 


3.6 x 10~ 12 above 250 GeV [46] 


Draco 


VERITAS 


20 


2.4 x 10~ 12 above 200 GeV [44] 


Ursa Minor 


VERITAS 


20 


2.4 x 10~ 12 above 200 GeV [44] 


Willman 1 


VERITAS 


15 


2.4 x 10~ 12 above 200 GeV [44] 


Willman 1 


MAGIC 


15.5 


10~ 12 above 100 GeV [47] 


Draco 


MAGIC 


7.8 


1.1 x 10~ n above 140 GeV [48] 


Draco 


STACEE 


10.2 


1.2 x 10- 7 (£/GeV)- 2 - 23 [42] 


Draco 


Whipple 


14.3 


E 2 (dF/dE) at 400 GeV < 5.10 x 10~ 12 erg[49] 


Ursa Minor 


Whipple 


17.2 


E 2 (dF/dE) at 400 GeV < 7.30 x 10" 12 erg[49] 



Table 3: Previous observations of dwarfs by ACTs and corresponding bounds on photon flux 



Figure 1 compares these bounds with the predictions of models A, B, C. The signal 



fluxes are computed using the astrophysical factors in Table 1. We assume the field of 
view corresponding to Ad = 0.15° in all cases. Experimental bounds from STACEE and 
Whipple, not shown in the figure, place weaker constraints due to smaller effective areas 
of the telescopes and lower signal fluxes from Draco and Ursa Minor. Figure 1 shows that 
the predictions are consistent with the observed null results within the uncertainties in the 
astrophysical factors. 



10 



VERITAS bound 
>200 GeV 



10 14 r Draco 



10 



Ursa Minor 



MAGIC on 
Willman 1 Willman 1 

>100GeV 



dwarfs and models 



MAGIC 
on Draco 
>140 GeV 



Figure 1: Comparison of experimental bounds with predictions from theory. The horizon- 
tal lines represent the bounds from Table 3. The three vertical bars for each search are 
the corresponding predictions of models A (left bar), B (center), and C (right), using the 
astrophysical factors from Table 1. 



3.4 Backgrounds 

Since the dwarf galaxies themselves are not expected to contain significant sources of hard 
gamma rays of astrophysical origin, the background can be effectively measured by looking 
at a region of the sky close to the dwarf (called the OFF region). Subtracting the OFF region 
flux from the flux in the ON region (which contains the dwarf) eliminates the background, 
up to statistical fluctuations, as long as it is isotropic on the relevant angular scales. If N 
background photons are expected in the ON region in a given energy bin, the subtraction 
method effectively reduces the background to \/N. In order to estimate this residual, which 
limits the sensitivity of dark matter searches, we use the standard extrapolations of the 
charged lepton, hadron and gamma-ray cosmic ray spectra [50]: 



rf 2 $ 



lep 



d 2 $ 



dEdVt 

had 



1.73 x 10" 



3.3 



dEdVt 



4.13 x lO"^, 



dEdtt 



3.6 x 10 



ad 



-10 



100 GeV\ 

E ) 
100 GeV\ 2 - 7 



cm 



E 



/100 GeV\ 



cm 



2.7 



V E 



) 



cm 



s- 1 GeV- 1 sr- 1 , 


(10) 


s- 1 GeV- 1 sr- 1 , 


(11) 


s- 1 GeV- 1 sr- 1 . 


(12) 




Figure 2: Dependence of instrument sensitivity on ehad, normalized to 1 at ehad = 1- The 
dashed, solid, and dot-dashed lines correspond to energy thresholds of 500, 200, and 100 
GeV respectively. 

Here ehad accounts for rejection of hadronic jets, normalized so that ehad = 1 for the Whipple 
telescope. We will conservatively set ehad = 1 throughout our analysis. Figure 2 shows 
how instrument sensitivity depends on ehad- Sensitivity improves with hadron rejection, 
but for sufficiently small values of ehad (below about 0.01) the lepton background starts to 
dominate, after which the exact value of ehad becomes irrelevant. The figure shows that, with 
improved hadron rejection, sensitivity can improve by up to a factor of 2 to 3 compared to 
the sensitivity with ehad = 1, depending on the threshold energy. 

Note that all backgrounds are assumed to be isotropic. We will ignore the gamma-ray 
contribution in our analysis since it is significantly smaller than the hadronic and leptonic 
backgrounds. It is worth keeping in mind that these formulas are simple power law extrapola- 
tions from observations at lower energies; we are only relying on these to make estimates, and 
actual observations will directly measure and eliminate this background through background 
subtraction, as mentioned earlier. 

In addition to astrophysical background, gamma rays due to dark matter annihilation in 
the Milky Way should be considered as background for dwarf searches as they contribute 
approximately equally to the flux in the ON and OFF regions. In the leptophilic dark mat- 
ter models we are interested in, these gamma rays come primarily from inverse Compton 
scattering (ICS) of starlight and CMB photons off of energetic positrons produced by WIMP 
annihilation, as well as the FSR process occurring in the Milky Way. We have calculated 
the ICS distribution using the semi-analytic formalism in [5] and [51]. An isothermal profile 
for the dark matter density, and the MED propagation model, where the half thickness of 
the galactic diffusion zone is 4 kpc, have been assumed. Both the ICS and FSR backgrounds 
depend on the line of sight from the Earth to the respective dwarf and need to be calcu- 
lated separately for each dwarf. In addition, these backgrounds depend on the dark matter 
annihilation channel. 

Table 4 lists the expected signal and background fluxes, integrated over energy above 200 



9 



Dwarf 


Draco 


Segue 1 


Segue 1 






(A0 = O.15°) 


(A0 = 0.10°) 


signals 








Signal Model A 


2.6 x KT 13 


9.8 x 10~ 12 


7.3 x 10~ 12 


Signal Model B 


9.9 x KT 14 


3.7 x KT 12 


2.7 x 10~ 12 


Signal Model C 


6.5 x KT 14 


2.4 x KT 12 


1.8 x 10~ 12 


backgrounds 








non-WIMP bg 


2.o x kt 11 


2.0 x 10~ n 


8.8 x 10~ 12 


ICS Model A 


3.3 x KT 16 


1.5 x 10~ 16 


6.8 x 10~ 17 


ICS Model B 


4.8 x KT 16 


2.2 x 10~ 16 


9.9 x 10~ 17 


ICS Model C 


2.5 x KT 16 


1.2 x 10~ 16 


5.2 x 10~ 17 


FSR Model A 


2.3 x KT 15 


1.3 x 10~ 15 


5.9 x 10~ 16 


FSR Model B 


8.5 x KT 16 


4.9 x 10~ 16 


2.2 x 10~ 16 


FSR Model C 


5.6 x KT 16 


3.2 x 10~ 16 


1.4 x 10~ 16 



Table 4: Signal and background fluxes (in cm s ) above 200 GeV for Draco and Segue 1. 
Dark matter-induced ICS and Milky Way FSR backgrounds are subdominant to astrophysical 
backgrounds and can be ignored. 



GeV, for the dwarfs Draco and Segue 1 for each annihilation model. While Draco has the 
lowest central value for its astrophysical factor among the dwarfs in Table 1 and is closest 
in direction to the galactic center, Segue 1 is the exact opposite; hence these provide two 
extreme cases. The signal fluxes were calculated using the central values of the astrophysical 
L factors listed in Table 1; in the final column of Table 4, we used the astrophysical factor 
for Segue 1 for an angular resolution of 0.10°, the appropriate parameter for MAGIC-II [16]: 



log in (L x GeV~- 



ilOV 



cm 



20.04 ±1.40. 



(13) 



The WIMP backgrounds are found to be negligible compared to the astrophysical back- 
grounds for all sources considered here, and will be ignored for the rest of the paper. Depend- 
ing on the source, the telescope, and the dark matter model, the signal/background ratios 
vary between about 1 and 0.005. For a higher threshold, say 500 GeV, the signal/background 
ratio will be higher since the background falls off faster than the signal at higher energies. 
As explained above, the background is measured directly from data, allowing for convincing 
discovery if sufficient statistics can be accumulated. 

4 Detection Prospects 

In this section, we will estimate detection prospects of current and near-future ACTs for FSR 
photons within the three dark matter models at hand. We will assume that the background 



10 



subtraction method is used, and ignore systematic errors. For positive detection, we impose 
the requirements that the significance of the excess in the ON region relative to the OFF 
region exceed 3a (or 5a), and the number of photons in the ON region exceed 25: 



Significance 



$ 7 A 



off 1 



(A eS ty/ 2 {d<s> hg /dQ x Any/ 2 



> 3 (or 5), 



Number of signal events = $ 7 A eff t > 25, 



(14) 



(15) 



where A er r is the effective area of the instrument and t is the observation time. For these 
estimates the angular sizes of the ON and OFF regions were assumed to be equal, and set to 
the minimum value consistent with the angular resolution. 7 It is found that for most ACT 
observations of dwarf galaxies, measurement is limited by backgrounds rather than statistics. 
To improve sensitivity, it is useful to choose an energy threshold that maximizes the 
ratio $ S ig na i / J <&bg- In Figure 3 we plot this ratio as a function of threshold energy for the 
three annihilation channels. The optimum energy threshold lies between 200 and 700 GeV 
depending on the model. For our estimates we will use two common values, 200 GeV and 
500 GeV, for all three models to allow direct comparison between models. 



1.2 



1.0 



•t 0.8 



0.6 




200 



400 600 800 1000 

Threshold Energy 



Figure 3: ^ S i gna i/ J&bg as a function of threshold energy above which the flux is integrated 
for the three annihilation channels of interest, normalized to 1 at 200 GeV. 

Figure 4 shows the length of observation time needed for a 3a detection of the FSR signal 
from each dwarf for each annihilation model, for two energy thresholds: 200 and 500 GeV. 
Likewise, Figure 5 shows the approximate minimum integrated flux above 200 GeV that can 
be detected at the 3a level in 50 hours of observation for a few ACTs, and the flux predicted 
by each model for various dwarfs. These sensitivities are also listed in Table 5, along with 
the corresponding number of events. 



7 Best sensitivity is achieved by minimizing the size of the ON region to improve S/VB, as long as the 
ON region is large enough to capture most of the FSR photons from the dwarf. Increasing the size of the 
OFF region is not desirable, since angular variations in the background may become an issue, and in any 
case the residual background is limited by the ON region size. 
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Figure 4: Observation times needed for 3a detection with MAGIC parameters. The six 
points for each dwarf correspond to (from left to right) models A, B and C with a 200 GeV 
energy threshold, and models A, B and C with a 500 GeV threshold. Vertical bars correspond 
to the uncertainties in the astrophysical factors of the dwarfs from Table 1, with the circles 
corresponding to the central values from the table. The solid horizontal line denotes 50 hours 
of observation time. The dashed line, at 500 hours, is equivalent to 50 hours of observation 
with an order of magnitude increase in the effective area, i.e. CTA parameters. 
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Figure 5: Integrated fluxes above 200 GeV. The dot-dashed, solid, and dashed lines cor- 
respond to approximate sensitivities of VERITAS, MAGIC, and CTA respectively, for 50 
hours of observation time (see Table 5). 

It is clear from the plots that the uncertainties in the astrophysical factors, and therefore 
observation times and sensitivities required for a positive signal, span several orders of mag- 
nitude; as a result, no models are ruled out, and no dwarf is guaranteed to give an observable 
signal. Draco has the smallest uncertainty in its astrophysical factor, but is also the least 
likely to give an observable signal. Segue 1 provides excellent prospects for detection, but the 
uncertainty on its astrophysical factor is huge. In summary, these estimates only allow us to 
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ACT 


Flux Sensitivity 
(> 200 GeV) 


Number of events 
(> 200 GeV) 


Flux Sensitivity 
(> 500 GeV) 


Number of events 
(> 500 GeV) 


VERITAS 

MAGIC-II 

CTA 


1.2 x 10-i 2 
9.9 x 10-13 
3.14 x 1Q- 13 


150 
180 
565 


5.1 x 10-13 

4.3 x 10-13 

1.4 x 10-13 


65 

77 
245 



Table 5: Minimum signal flux (in cm"^" 1 ), and the corresponding number of signal events, 
needed for 3<t detection of an excess above the astrophysical background with 50 hrs of 
observation time. 



conclude that, for current and future ACTs, detection of FSR from dark matter annihilation 
from the above dwarfs is likely, but not guaranteed. Further astronomical observations of 
the dwarfs should reduce the uncertainty in the L factors, allowing more precise predictions 
to be made in the future. 

Throughout this section, we conservatively set e^ad = 1, corresponding to the hadron 
rejection capabilities of the Whipple telescope. Fig. 2 shows that the sensitivity of the 
searches can only be enhanced by factors of 2-3 with improved hadronic rejection, before 
the lepton background becomes dominant. Thus, while certainly desirable, hadronic rejection 
improvements will not affect the qualitative conclusions of our analysis. 

5 Model Identification: Some Case Studies 

In this section we investigate the prospects of identifying the dark matter model based on an 
observed FSR gamma ray signal. As explained above, by "model" we really mean the final 
state of WIMP annihilation; to leading approximation used throughout this study, the FSR 
signal is completely insensitive to details of the microscopic model giving rise to WIMPs and 
their annihilation. For concreteness, we will focus on the "models" A, B and C, defined in 
Section 2. 

Our approach is as follows: We assume that one of the models, with particular parameter 
values, is realized in nature. (The parameters are (crv), m x , and, in model C, m^. Note that 
m,p is not an independent parameter for Model B, where it can be absorbed into (crv) since 
there is no additional contribution to FSR from the decay of the final products; see Eq. (4).) 
We generate a set of random "data points" distributed in energy in accord with the theoretical 
predictions of this model, incorporating the energy resolution of the instrument according 
to Eq. (9). The total number of data points corresponds to the prediction of the model for 
a particular source and telescope parameters. These data points are then binned, choosing 
the bin width to be approximately double the energy resolution of the ACT at hand. The 
WIMP mass will in practice be unknown, but the final bin can always be made large enough 
that E 1 = m x can be assumed to fall in this bin. 

For each data set, two background samples, corresponding to the ON and OFF regions, 
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are then generated. The number of background events in each bin is found as the difference 
between the event counts in the ON and OFF samples in that bin. A fit to the binned 
data points is performed with all three models 8 , varying the parameters in each model to 
minimize x 2 jd.o.j .. The three optimal x 2 /d.o.f. values obtained in this way are compared 
to each other, and the model with the smallest value is declared the "best fit" to the data. 
We perform this procedure for 100 randomly generated data sets, and collect information on 
the number of data sets for which each of the models (A, B and C) provides the best fit, as 
well as on the best-fit values of parameters. For most cases, we verified the convergence of 
model discrimination, as well as the best fit mass and cross section, with a larger number of 
data sets (about 300). 

We first perform the above analysis assuming that the model A is true, and then repeat 
it assuming that the "true" model is B and C. When generating the data set, we fix the 
model parameters to m x = 3 TeV, (av) = 1000 pb, and, where relevant, m^ = 1 GeV. We 
perform the analysis for two sets of ACT parameter values: A e g = 10 9 cm 2 and 15% energy 
resolution, as for MAGIC-II (which we label "current" below, since these are representative 
of the reach of current instruments); and A e g = 10 10 cm 2 and an energy resolution of 10%, 
as for the CTA (which we label "future"). We use A6 = 0.10° for both ON and OFF regions 
in all cases. 

The results of the analysis are summarized in Tables 6-10. Each table contains the 
distribution of "best fit" models, best fit WIMP mass, cross section, and x 2 /d.o.f.. (Fits 
to rriff, have very large error bars in almost all cases and therefore are not listed.) The first 
entry, the distribution of "best fit" models, contains three numbers, corresponding to the 
number of data sets for which models A, B, and C respectively gave the best fit. For the 
other parameters we also list the statistical error bar reflecting the variation of the best-fit 
values and the x 2 /d.o.f . of the best fit among the 100 data sets. Discussion of success rates 
of model identification is based on the assumption that all three models are a priori equally 
likely. In practice, information from independent measurements {e.g. positron spectra) can 
be used to place priors, enhancing the model discriminating power of the FSR signal. 

Case 1: Observation of Segue 1 

As our first, "benchmark" case study, we consider the observation of the most promising 
candidate Segue 1, assuming the astrophysical L factor at the current central value from 
Eq. (13) and an observation period of 50 hours. We consider a data sample with an energy 
threshold of 500 GeV, distributed in 6 (8) energy bins for current (future) telescope param- 
eters. The results of this analysis are summarized in Table 6, and a sample fit to one of the 
100 data sets is shown in the left panel of Fig. 6. Note that the cross section errors listed 
in Table 6 and below do not include the uncertainty of the astrophysical factor L. In other 
words, the errors reflect the relative accuracy of the measured normalization of the signal 
flux, which is proportional to the product o x L. 



8 We use analytic predictions of the models in these fits, ignoring energy resolution of the telescopes. 
Given the chosen bin widths, energy resolution effects are unlikely to significantly affect the fits in our study. 
In a more complete analysis, the fit should be performed using Monte Carlo simulations of the signal in each 
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"True" model 


model A,B,C 


Best fit WIMP 


Best fit 


X 2 /d.o.f. 




as best fit 


mass (GeV) 


cross-section (pb) 




current 










model A 


too, 0, 


2978 ±112 


1013 ±34 


1.16 ±0.75 


model B 


16, 44, 40 


3083 ±412 


1054 ± 129 


1.20 ±0.95 


model C 


17, 13, 70 


3073 ± 449 


1199 ± 506 


1.51 ±1.11 


future 










model A 


100, 0, 


3006 ± 28 


999 ±9 


1.11 ±0.55 


model B 


0, 68, 32 


3097 ± 144 


996 ± 38 


1.37 ±0.89 


model C 


0, 11, 89 


3085 ± 389 


983 ± 269 


1.42 ±1.05 



Table 6: Case 1 fit results. 

These fits provide several noteworthy results, which also apply to all variations of the 
analysis considered below (Cases 2-5): 

• The 2/i channel is very clearly identified. The two 4-body final state channels are 
not easily distinguishable from each other because their spectra are very similar, see 
Figure 6. The key distinguishing feature of the 2/z channel is the sharp edge at the end 
of the spectrum, common to all two-body fermionic channels [17]. In contrast, 4-body 
channels have a softer spectrum, gradually falling off to zero at high energies. 

• Between the 4-body channels, the 4// channel is more likely to give a better fit due to 
the extra degree of freedom in m^. For this reason, it is often preferred even if the true 
model is 4e. 

• Measured values of the WIMP mass and annihilation cross section are in excellent 
agreement with the "true" values. The mass is better reconstructed with future tele- 
scope parameters because of superior energy resolution, as expected. The accuracy of 
the mass and cross section determination is particularly impressive (of order 1%) in 



the 2/x channel, again presumably due to the sharp edge feature at E 1 



in 



x- 



• Future telescope parameters provide a clear improvement over current parameters in 
terms of correct identification of the model. This is mainly due to the increase in 
effective area, which results in a greater number of events and a better significance. 

Overall, model identification had a success rate of 71% for current telescope parameters 
and 86% for future parameters in this case. In the cases that follow, we examine the effects 
of changing some of the assumed parameters. 

Case 2: A lower threshold 



model, which would include energy resolution and other instrumental effects. 
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Figure 6: Left: A sample fit. In this particular example, the data is generated according to 
the model A, which is also correctly identified as the best-fit model. Right: Energy spectra 
from theory. The thick, dashed, and thin lines correspond to models A, B, and C respectively. 
The WIMP mass is assumed to be 3 TeV in all models. The flux has been normalized to 1 
at 500 GeV for all three curves. 

To estimate the effect of lowering the energy threshold on model-discriminating power, 
we repeated the benchmark analysis for a threshold energy of 200 GeV, with data distributed 
into 9 energy bins (13 for future parameters). All other parameters are the same as in Case 1. 
The results in Table 7 show that there is overall improvement in model identification, but the 
X 2 /d.o.f. values get slightly worse due to higher background flux at lower energies. Model 
identification success rates are 79% and 93% for current and future telescope parameters 
respectively. Based on this, there might be some reason to prefer a lower energy threshold, 
but performing multiple fits with multiple energy threshold cuts on the same signal might 
be a better approach to identifying the correct model. 



"True" model 


model A,B,C 
as best fit 


Best fit WIMP 

mass (GeV) 


Best fit 
cross-section (pb) 


X 2 /d.o.f. 


current 
model A 
model B 
model C 


99, 1, 
12, 59, 29 
5, 16, 79 


2987 ±107 
3125 ±372 
3137 ±436 


1009 ±26 
1024 ± 72 
1062 ±367 


1.28 ±0.69 
1.38 ±0.83 
1.53 ±0.76 


future 
model A 
model B 
model C 


100, 0, 
0, 83, 17 
0, 3, 97 


3003 ± 24 
3074 ± 134 
3155 ±347 


1000 ± 7 

997 ± 28 

1018 ±183 


1.44 ± 0.62 

1.52 ±0.67 

1.53 ±0.71 



Table 7: Case 2 fit results. 



Case 3: Improved hadron rejection 



16 



Large atmospheric background is a significant hurdle for dark matter detection from 
dwarf galaxies using ACTs, as well as for model discrimination. In the benchmark analysis 
we conservatively set e^ad = 1, corresponding to the hadron rejection capabilities of the 
Whipple telescope. Although we do not have exact values, newer instruments, both current 
and future, are likely to possess significantly improved hadron rejection capabilities. For this 
case we consider e^ad — 0.01. According to Figure 2, this should improve sensitivity by a 
factor of 3, and further reducing e^ does not have much of an effect since the background 
becomes dominated by the lepton flux. The remaining parameters are the same as in Case 
1. The results, listed in Table 8, show some improvement in model identification in all cases 
. Overall, hadron rejection capabilities of a telescope do not seem to be a limiting factor for 
model discrimination. 



"True" model 


model A,B,C 

as best fit 


Best fit WIMP 
mass (GeV) 


Best fit 
cross-section (pb) 


X 2 /d.o.f. 


current 
model A 
model B 
model C 


100, 0, 
5, 53, 42 
0, 6, 94 


2999 ± 97 
3082 ± 269 
2951 ± 366 


1005 ±28 
1014 ±68 
575 ± 366 


1.06 ±0.74 
0.90 ±0.87 
1.15 ±0.78 


future 
model A 
model B 
model C 


100, 0, 
0, 74, 26 
0, 0, 100 


3007 ± 26 
3080 ± 83 
3021 ± 238 


999 ±7 

987 ± 20 

1011 ±184 


1.04 ±0.58 
1.06 ±0.59 
2.32 ±2.68 



Table 8: Case 3 fit results. 

Case 4: Lighter dark matter 

The benchmark analysis assumed a WIMP mass of 3 TeV, which is at the higher end of 
masses that give good fits to PAMELA, Fermi and HESS data. While a lower mass results in 
fewer energy bins and/or comparatively larger background, it also leads to a larger signal due 
to the m~ 3 dependence of the FSR signal in Eqs. (1), (4). To estimate the effect on model 
discrimination, we repeated the benchmark analysis with m x = 1 TeV, an energy threshold 
of 200 GeV, and 5 energy bins (7 for future parameters). The results, listed in Table 9, 
are very similar to the results from case 2, which uses the same energy threshold; the only 
notable difference is the improved ability to clearly distinguish 2-body channels from 4-body 
channels with current telescope parameters. Thus, model discriminating capabilities of the 
ACTs are generally quite robust across the range of WIMP masses preferred by current data 
within leptophilic models. 

Case 5: 3<r or 5a detection 

In all studies above, we have used the central value of the astrophysical factor L of Segue 
1, which results in rather high fluxes that can be easily detected by MAGIC and CTA, see 
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"True" model 


model A,B,C 


Best fit WIMP 


Best fit 


X 2 /d.o.f. 




as best fit 


mass (GeV) 


cross-section (pb) 




current 










model A 


100, 0, 


997 ±16 


1004 ±16 


1.02 ±0.77 


model B 


0, 51, 49 


1057 ±62 


980 ± 42 


1.05 ±0.88 


model C 


0, 4, 96 


1064 ± 123 


997 ± 354 


1.31 ±1.65 


future 










model A 


100, 0, 


1000 ±4 


1001 ±4 


1.26 ±0.81 


model B 


0, 62, 38 


1021 ± 16 


985 ± 14 


1.28 ±0.74 


model C 


0, 0, 100 


1040 ± 62 


1031 ±130 


1.21 ±0.77 



Table 9: Case 4 fit results. 

Fig. 5. Due to large uncertainties in the astrophysical factors, in reality fluxes may be much 
lower. For the lowest possible values of L, the ACTs may not be able to even detect the FSR 
signal from Segue 1; of course, no model discrimination would be possible in this case. It is 
interesting, however, to examine how much can be learned about the underlying model if the 
FSR flux is just borderline observable. We performed this analysis, assuming the weakest 
possible signals that can be detected at 3a or 5a significance. All other parameters are as in 
case 1. Table 10 shows the fit results for current telescope parameters. (Not surprisingly, fit 
results for future telescope parameters, for the same observation significance levels, are very 
similar, and therefore not shown.) 



"True" model 


model A,B,C 

as best fit 


Best fit WIMP 
mass (GeV) 


Best fit 
cross-section (pb) 


X 2 /d.o.f. 


3<7 detection 
model A 
model B 
model C 


71, 12, 17 
2, 62, 36 
0, 49, 51 


3092 ±531 
3041 ± 524 
3051 ± 503 


1235 ± 265 

1236 ± 275 
1238 ± 650 


1.17 ±0.71 
1.74 ±1.62 
1.26 ±0.77 


5<7 detection 
model A 
model B 
model C 


83, 9, 8 

2, 64, 34 

3, 17, 80 


3146 ± 433 
3117 ±471 
3065 ± 489 


1087 ±182 
1148 ±238 
1005 ± 629 


1.32 ±0.76 
1.51 ±0.97 
1.44 ±0.95 



Table 10: Fit results for the weakest possible signal detectable at 3(7 and 5(7 levels by the 
current ACTs. 

The fit results show that the overall success rate for correctly identifying the model in 
these scenarios is comparable to the benchmark case 1: 61% and 75% for 3<r and 5a detection 
respectively. The success rates of correctly identifying the signal as a 2-body or 4-body final 
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state remain excellent: 90% and 93% respectively. Best fit values obtained for the WIMP 
mass and annihilation cross section remain fairly accurate. We stress again that these results 
are for the weakest signals that can be detected at the given significance, and any stronger 
signal will give better fits and model identification. 

6 Conclusions 

The following is a summary of our main findings: 

• We reviewed the prospects of indirect detection of dark matter via observation of final 
state radiation (FSR) gamma rays from dwarf galaxies using current and near-future 
atmospheric Cherenkov telescopes. We found that, assuming dark matter annihilation 
into leptonic final states, with parameters motivated by the recent PAMELA, Fermi 
and HESS anomalies, a detection of such a signal is very likely. Unfortunately, lack 
of precise knowledge of the distribution of dark matter in the dwarfs makes the signal 
flux predictions highly uncertain, so that null results of these searches would not con- 
clusively rule out these models. (In other words, positive detection is not guaranteed.) 

• If a signal is observed, fits to the photon spectrum can be used to discriminate between 
dark matter models. Prospects are especially encouraging for the most promising can- 
didate dwarf Segue 1. With the current central value of its astrophysical factor L, our 
simulations show a 71% (86%) success rate of correctly identifying the dark matter 
model with current (future) telescope parameters. The success rate for determining 
whether the annihilation is into a 4-body final state or a 2-body final state - a dis- 
tinction that can be tremendously helpful for restricting viable models of dark matter 
- is close to perfect. Distinguishing between 4e and 4/j, channels is more challenging, 
because their energy spectra are very similar, but future telescopes are expected to 
provide a significant improvement in this regard. 

• Even if a signal is barely detected at a 5a level, we estimate that there is at least a 
75% chance of correctly identifying the annihilation channel among the three options 
discussed in this paper. In particular, if annihilation is into two muons, there is about 
an 83% chance of correctly identifying it. 

• Fits to photon energy spectrum can also be used to determine the values of the WIMP 
mass and annihilation cross section (or, more precisely, a x L). For a signal detected 
at the 5<T level, these parameters can be determined with roughly 15% accuracy. Note 
that no prior knowledge of the model is assumed in these fits. 

• The success rate for proper model identification is fairly robust with respect to changes 
in energy threshold, WIMP mass, energy resolution, and hadron rejection capabilities 
of the telescope, and seems to depend mainly on the S/y/~B ratio. 
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In summary, if PAMELA, Fermi and HESS anomalies have their origin in leptophilic dark 
matter annihilation, current and near-future ACTs have an excellent chance of observing the 
accompanying FSR gamma ray signal from dwarf galaxies. Once a signal is observed, the 
measured gamma ray spectrum can be used to identify the correct annihilation channel, 
paving the way to a better understanding of the microscopic nature of dark matter. Our 
study shows that this approach is very promising. It should be kept in mind, of course, 
that the phenomenological analysis presented here does not take into account potentially 
important instrumental effects specific to each telescope, such as instrumental backgrounds, 
systematic errors on photon energy measurement, energy dependence of the effective area, 
etc. We hope that the positive results of our study will encourage experimental collaborations 
to perform a more careful examination of their instruments' capabilities in this regard, and 
to continue searches for dark matter annihilation signals from dwarf galaxies. 

Note Added 

Since the submission of this manuscript, further analysis has driven the astrophysical 
factor estimates for Segue 1 to lower values [52]; however, we have verified that the general 
conclusions of our analysis still hold. In particular, the results for Case 5 in Section 5, 
which depict worst case scenarios with positive detection, are independent of the value of 
the astrophysical factor of Segue 1. 
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